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We discuss a possibility to directly reconstruct the CMB polarization field at the last scattering 
surface by accounting for modifications imposed by the gravitational lensing effect. The suggested 
method requires a tracer field of the large scale structure lensing potentials that deflected propagating 
CMB photons from the last scattering surface. This required information can come from a variety 
of observations on the large scale structure matter distribution, including convergence reconstructed 
from lensing shear studies involving galaxy shapes. In the case of so-called curl, or B,-modes of 
CMB polarization, the reconstruction allows one to identify the distinct signature of inflationary 
gravitational waves. 
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The advent of high sensitivity cosmic microwave background (CMB) experiments will soon allow detailed observa- 
tional studies on the CMB polarization field. In addition to a confirmation of the basic picture on how anisotropies 
were generated during the recombination era Q, CMB polarization observations, at large angular scales, may allow a 
measurement of the reionization rcdshift based on excess power created during the rescattering of CMB quadrupole 
at the reionization surface || . An additional use of polarization observations is the potential detection of inflationary 
gravitational waves based on its contribution to the so-called curl or B-modes of polarization |3j . 
t-H ■ In addition to primordial contributions involving density perturbations, which contribute solely to E-modes, and 
tensor modes or gravitational waves, the polarization field one observes today also contains secondary contributions 
associated with the large scale structure. This is similar to well known secondary contributions that dominate the 
^{■^ small scale temperature anisotropy signal Q. In the case of B-modes, the arcminute scale polarization is dominated 
iy~j • by a fractional conversion of the dominant E-mode contribution via the gravitational lensing angular deflection of 
CMB photons f§. In Ref. §, we discussed an approach to separate the gravitational-wave signature in B-modes 
from the dominant, and contaminant, weak lensing contribution at arcminute scales. We also investigated limits 
on gravitational wave B-mode detection after a model-independent lensing substraction. This calculation followed 
■ previous discussions in the literature on how CMB data can be used for an extraction of statistics involving the 
Q-f gravitational lensing effect, such as shear or deflection angles [^[-^). 

In general, the reconstruction of lensing deflections from CMB data require a priori information on the primordial 
CMB contribution at the last scattering surface. This knowledge cannot easily be obtained from data alone since 
the CMB contribution measured today involve contributions resulting from secondary effects and any modifications 
during the transit to us from the last scattering surface. On the other hand, we can consider an alternative approach 
to an analogous problem involving a reconstruction of the primordial CMB contribution, at the last scattering surface, 
by accounting for modifications due to large scale structure. This is, effectively, an inverse problem from what was 
discussed in [|||] and require prior information related to the lensing effect on CMB, such as the mass distribution 
of large scale structure in which photons propagated. In the present study, we briefly discuss this possibility using a 
tracer field of the large scale structure potentials that deflected CMB photons via the gravitational lensing effect. 

The discussion presented here follows the recent work of Ref. § where we discussed the lensing reconstruction from 
CMB data given information related to primordial CMB anisotropy and polarization fields. We refer the reader to 
the above paper for basic details of lensing and other important ingredients related to the calculation. To illustrate 
our results, we assume the currently favored ACDM cosmology. 



II. CALCULATION 



The lensing effect on CMB data can be described simply as a transfer of power via changes to photon propagation 
directions on the sky. In Fourier space, we can consider the E- and B-mode decomposition introduced in Ref. || such 
that ±Y(l) = E(V) ± iB{\). We denote the primordial contribution at the last scattering surface with a Y, while the 
polarization field observed today is denoted with Y. We can consider the lensing modification following Ref. and 
write 

± y(l) = ± f(l)-| -^^(lJe^H-^L&h), (1) 
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where 
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to the second order in the projected lensing potential <fi given by [|5]Jl 
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Here, r is the conformal time or comoving radial distance, d A is the analogous comoving angular diameter distance, $ 
is the gravitational potential. The latter is related to fluctuations in the density field via the Poisson equation. The 
deflection angle associated with lensing is the gradient of the projected potential, a(n) = V0(n). Using deflections 
to the first order, we write 



E(l) = E(l) 
B(l) = B(l) 
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In order to extract the primordial contribution, unaffected by lensing, consider the combination of E and B fields 
with a tracer field of the deflecting potential, which we label here as X. This tracer map has the property that it does 
not correlate strongly with either E or B, the quantities at the last scattering. Thus, only correlation expected is the 
one resulting due to lensing, and to a lesser extent from other secondary contributions to polarization. Since latter 
contributions involving scattering effects associated with large scale structure is significantly smaller jll]], we will 
ignore them in this discussion. Additionally, secondary scattering effects associated with polarization can be removed 
from data based on their spectral dependences. This is similar to the approach suggested in Ref. (T^| to remove the 
Sunyaev-Zel'dovich effect (SZ; from dominant temperature anisotropies. We consider a map created by taking 
the product of the polarization field and the tracer field, e.g., EX . In Fourier space, we can write this product as a 
convolution of the Fourier moments of the two fields: 



{EX){\) = \ J ^±[E(h)X(l-h)+X(h)E(l-h)]W(l,h 
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where W is a Fourier-space based filter that we will design to maximize the reconstruction of the primordial polarization 
information. Similarly, we can also construct another quadratic statistic involving B and X-fields. Using equation (|^) 
in (H), we can simplify to obtain 
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while, analogous contribution to the quadratic statistic involving B and X fields, (BX)(l), is 



(BX)(l) = Ie(1) 
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As written, these quadratic combinations are proportional to the primordial polarization contributions at the last 
scattering surface, E(l) and B(l). This is essentially the basis for the suggested reconstruction. In suggested 
quadratic statistics involve combinations of EE and EB such that these lead to terms which are proportional to 
<fi times integrals involving Cf E and Cf B \ Thus, with some knowledge on these contributions, it was possible to 
reconstruct lensing deflection potentials. Our approach, though similar, involves the inverse problem of constructing 
E and B given information related to how <fi behaves. In above, Cf^ is the cross-correlation between tracer X field 
and the lensing potentials <f> defined such that 
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(X(lW)} = (27T) 2 fc(l + l')Cf* . 



We can write this cross angular power spectrum as 
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and similarly for lf(k). Here, W l (k, r) is the window function of the X or (f> in radial coordinates, ji's are the spherical 
Bessel functions and Pgg(k) is the power spectrum of density fluctuations (see, Ref. || for details). 

To extract power spectra associated with the primordial polarization, Cf E and C BB , we can either use power 
spectra constructed from correlating quadratic statistics themselves, e.g., {(EX)(Y)(EX)(Y)), or the quadratic statistic 
correlated with the polarization field, e.g., ((EX)(l)E(Y)) . The latter approach involves an extraction of the power 
spectrum from the bispectrum formed by the fields traced by the quadratic statistic and either E or B; this is similar 
to what was suggested in Ref. to extract the lensing-SZ correlation from squared temperature-temperature power 
spectrum. The method based on ((EX)(V)(EX)(Y)) extracts the power spectrum information via the trispectrum 
formed by lensing and tracer field correlation associated with quadratic statistics. This is analogous to techniques 
in Refs. ||[| to extract the lensing information using CMB data. Here, we will use the approach based on the 
bispectrum since it has less noise contributions; the approach based on the trispectrum has additional non-Gaussian 
noise contributions which can complicate the extraction while lowering the associated signal-to-noise ratio || . 

We first consider the power spectrum formed by ((EX)(Y)E(Y)) . Following our standard definition, we define the 
associated power spectrum in this case as 



((EX)(l)E(Y)) = (2n) 2 S D (l + Y)C E X ~ E . 

Similarly, we define the power spectra involved with other combinations of E and B with (EX) and (BX). 
We can write the angular power spectrum associated with the ((EX)(Y)E(Y)) correlation as 



(10) 



c, 



EX~E 



-C 



EE 



d 2 h 



C\ 



(1 - li) • I cos2(</?i -ip\i) + h ■ 1 cos2(y?i - ip\i-i t \) 



(11) 



Note that there is no contribution coming from the B term associated with (EX) (I) since there is no correlation 
between E(Y) and B(Y) due to parity considerations. Also, there is no resulting contribution associated with the first 
order lensing term of E as lensing deflections, (f>, should not correlate with the primordial polarization field. Similarly, 
the contribution to ((EX)(Y)B(Y)) results in a power spectrum given by 
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Since (BX) provides an addition statistic to correlated with polarization maps, we can also consider the cross- 
correlation of (BX)(Y) with E(Y) and 5(1). We can write the two contributions as 
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respectively. As written, C, E and Cf x E are proportional to Cf E while C EX B and Cf x B are proportional 
to C BB . This is the basis of the present approach which is leading to an extraction of these primordial polarization 
power spectra uncontaminated by the lensing contribution. With an appropriate normalization and a description 
for the filter, W(l, li) that optimizes the polarization extraction, one can effectively extract C BB and Cf E from the 
higher order statistics involving polarization maps today and a map of the large scale structure. 
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FIG. 1. Left: The CMB polarization B-modes power spectra. We show the error bars associated with reconstruction of the 
primordial power spectrum using no-noise E and B maps and a tracer of the large scale structure involving lensing potentials. 
For comparison, the dotted-line shows the contaminant contribution arising from the transfer of power from primordial E to 
B-modes via lensing deflections. Right: The cumulative signal-to-noise ratios for the reconstruction of the primordial power 
spectra. The long-dashed line is the signal-to-noise ratio associated with the detection of B-mode power spectrum, under the 
assumption that contaminant lensing contribution can be treated as a source of noise. 



The signal-to-noise ratio for the extraction of the primordial power spectra of polarization, for example, in the case 
of Cf B extracted using {(BX)(l)B(V)}, is 



/sky(2Z + 1) 



J $frW&h) {C$* h \ (! - !i) ' 1 cos2 (^ " Wx) + Cjf * li ' 1 cos2(^ - p Ml |)} 



NP X C, 



(15) 



where / s k v is the fraction of sky covered and noise contributions given by N, BX where 
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is the noise associated with (BX) statistic. In above, C; 3 '*, and C^'*, are the total contributions to B-mode and 
X-fields respectively. We can write these contributions as C\' = C\ + C" + Cf where Cf is any noise contribution, 
for example, instrumental noise in the B-map, and Cf any secondary contribution, such as the dominant lensing 
contribution in the B-map. 

In order to maximize the signal-to-noise ratio for the extraction of primordial polarization power spectra, we follow 
the approach in Ref. Pj and consider a description for the filter. The maximum signal-to-noise ratio is obtained when 
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We can define similar filters for the extraction of E from (BX) or (EX) estimators as well. In analogous to above, 
these filters involve the noise contributions in the (BX) and (EX) maps and the behavior of the function that 
integrates over the filter in each of the integrals in equations (|ll]) to (|l4|). Since each estimate of Cf E and Cf B 
comes from two independent statistics involving (EX) and (BX), we add the individual signal-to-noise estimates of 
the two. This is equivalent to weighing individual noise contributions inversely. 
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III. DISCUSSION 



We summarize our main results in figure 1 for the extraction of primordial polarization power spectra based on 
information related to lensing modification to CMB anisotropies. For the reconstruction, we assume a gravitational- 
wave contribution to B-modes with a tensor-to-scalar ratio of 0.25 and assume no instrumental-noise contributions 
to polarization observations today. When calculating noise, we include additional secondary noise contribution from 
lensing to E and B maps. For simplicity, we also take a tracer field, X, that corresponds to lensing potentials 
themselves. Thus, signal-to-noise ratios shown in figure 1 should be considered as optimistic values. Any noise 
contribution to polarization observations and a tracer field which is less correlated with <f> will lead to a degradation 
in the signal-to-noise ratios of primordial polarization reconstruction from the ones shown here. 

With X = 0, the assumption here is that one can use a properly cleaned CMB temperature map to reconstruct 
lensing deflections following the approach in MM and apply such a construction to polarization observations to reduce 
contamination from lensing to gravitational-wave detection. We include a noise-contribution to the X field, when 

X t n n ■ ■ ■ 

calculating C ; : , following the noise calculation for </> extracted from a no instrumental-noise CMB temperature map 
in Ref. |^,^). Additionally possibilities for X include a frequency cleaned SZ map, a map of the convergence from 
large scale structure lensing observations, the galaxy distribution from surveys such as Sloan but imaging out to a 
higher redshift; Since contributions to lensing in CMB comes at redshifts of a few and more, in general, a tracer field 
out to redshifts of ^ 3 to 5 will be required to appropriately apply this method. 

As suggested, the primordial polarization construction is more important for B-modes, given that they contain the 
distinct signature of gravitational waves, than E-modes. For the extraction of the B-mode power spectrum, we estimate 
a cumulative signal-to-noise ratio of ~ 16. The cumulative signal-to-noise for the detection of this contribution, under 
the assumption that the contaminant lensing contribution can be treated as a known source of noise, is of order ~ 
40. The suggested reconstruction here leads to a roughly a factor of ~ 2.5 decrease in the associated signal-to- noise 
ratio. The reconstruction, however, has the advantage that one does not need to make any assumptions on the lensing 
contribution to the power spectrum of B-modes. Relaxing our assumption of a tensor-to-scalar ratio of 0.25, we find 
that the suggested extraction can be used to separate primordial B-modes with at least a signal-to-noise ratio of 1 
when the tensor-to-scalar ratio is of order 0.01. Though we have not separated individual signal-to-noise values, the 
reconstruction of B-modes has less noise with the combination of (BX) statistic correlated with a map of B-modes 
instead of (EX) correlated with B-modes. This is because the dominant noise associated with E-modes in (EX) limit 
the reconstruction of .B-modes significantly. The signal-to-noise values for the reconstruction of .E-modes have similar 
values whether (EX) or (BX) is used. This reconstruction has a cumulative signal-to-noise in excess of 10 3 both due 
to its high contribution and the fact that E-mode power spectrum peak at arcminute scales compared to .B-mode 
power spectrum which peaks at degree scales. 

It is clear that the reconstruction analysis can be best done for the extraction of primordial -E-modes. The low signal- 
to-noise values for .B-modes, however, should not discourage from attempting to do analysis like the one suggested here 
in upcoming CMB data. The ultimate, and challenging, detection of gravitational wave signature, which, in principle, 
gives a handle on the inflationary energy scale requires various statistical techniques so that this contribution can 
be identified confusion free. Since lensing contribution dominates the B-mode signal, we have exploited its use as a 
possible way to extract primordial polarization field. Given the importance of detecting gravitational waves for future 
cosmological studies, statistics such as the one discussed are clearly warranted for further studies. 
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